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Abstract. We explore the viability of a bulk viscous matter-dominated Universe to 
explain the present accelerated expansion of the Universe. The model is composed 
by a pressureless fluid with bulk viscosity of the form ( = ( + (iH where Co an d 
Ci are constants and H is the Hubble parameter. The pressureless fluid characterizes 
both the baryon and dark matter components. We study the behavior of the Universe 
according to this model analyzing the scale factor as well as some curvature scalars 
and the matter density. On the other hand, we compute the best estimated values 
of Co and Ci using the type la Supernovae (SNe la) probe. We find that from all the 
possible scenarios for the Universe, the preferred one by the best estimated values of 
(C01C1) is that of an expanding Universe beginning with a Big-Bang, followed by a 
decelerated expansion at early times, and with a smooth transition in recent times to 
an accelerated expansion epoch that is going to continue forever. The predicted age 
of the Universe is a little smaller than the mean value of the observational constraint 
coming from the oldest globular clusters but it is still inside of the confidence interval 
of this constraint. A drawback of the model is the violation of the local second law 
of thermodynamics in redshifts z > 1. However, when we assume Ci = 0, the simple 
model C = Co evaluated at the best estimated value for Co satisfies the local second law 
of thermodynamics, the age of the Universe is in perfect agreement with the constraint 
of globular clusters, and it also has a Big-Bang, followed by a decelerated expansion 
with the smooth transition to an accelerated expansion epoch in late times, that is 
going to continue forever. 



PACS numbers: 95.36. +x, 98.80.-k, 98.80.Es 
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1. Introduction 

In recent years the type la supernovae (SNe la) observations have indicated a possibly 
late time accelerated expansion of the Universe [I]- pp. This discovery has been supported 
also by indirect observations of the expansion of the Universe like the cosmic microwave 
background (CMB) [8] and the large scale structure (LSS) [9] observations. It has been 
coined the name dark energy to refer to the unknown responsible of such acceleration 
that presumably is the ~ 70% of the total content of matter and energy in the Universe 

Since the time of this discovery in 1998 [1], a large amount of models have been 
proposed to explain this acceleration. The two most accepted dark energy models are 
that of a cosmological constant (assumed possibly to be the quantum vacuum energy) 
and a slowly varying rolling scalar field (quintessence models) |10]-|14j. 

Despite the cosmological constant model is preferred by the observations, this has 
several strong problems, among them, the huge discrepancy between its predicted and 
observed value (a difference of about 120 orders of magnitude) [15]- [18] and the problem 
of that we are living precisely in a time when the matter density in the Universe is of 
the same order than the dark energy density (the "cosmic coincidence problem") [llj, 
[191120]. 

On the other hand, in the context of inflation of the very early Universe, it has 
been known since long time ago that an imperfect fluid with bulk viscosity in cosmology 
can produce an acceleration in the expansion of the Universe without the need of a 
cosmological constant or some inflationary scalar field [21] [23] (although some authors 
do not agree with this conclusion [24J). So, extrapolating this theoretical idea used to 
induce an accelerating Universe without the need of unknown components it is possible 
to postulate that one candidate to explain the present acceleration can be a bulk viscous 
pressure [25]- [31]. However, this idea faces some problems too, for example, the need 
to have a satisfactory mechanism for the origin and composition of the bulk viscosity, 
nevertheless, there are already some authors working in this aspect [32j-[35j. 

From a thermodynamical point of view the bulk viscosity in a physical system is 
due to its deviations from the local thermodynamic equilibrium [31] (for a review in 
the theory of relativistic dissipative fluids see [36]). In a cosmological fluid, the bulk 
viscosity may arise when the fluid expands (or contracts) too fast so that the system 
does not have enough time to restore its local thermodynamic equilibrium and then it 
arises an effective pressure restoring the system to its thermal equilibrium. The bulk 
viscosity can be seen as a measurement of this effective pressure. When the fluid reaches 
again the thermal equilibrium then the bulk viscous pressure vanishes [31], [37]-|39j. 

So, in an accelerated expanding Universe, it may be natural to assume the possibility 
that the expansion process is actually a collection of states out of thermal equilibrium 
in a small fraction of time giving rise to the existence of a bulk viscosity. 

In the present work we study and test a matter-dominated cosmological model 
with bulk viscosity as an explanation for the accelerated expansion of the Universe. 
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The model is composed by a pressureless fluid (dust) with bulk viscosity of the form 
C = Co + QiH where Co an d Ci are constants to be determined by the observations and H 
is the Hubble parameter. The pressureless fluid characterizes both the baryon and dark 
matter components. The term "Co" takes into account the simplest parametrization for 
the bulk viscosity, a constant. And the term Ll (iH" characterizes the possibility of a 
bulk viscosity proportional to the expansion ratio of the Universe. 

The idea of this model is to drive the present acceleration using the trigger at recent 
times of the bulk viscous pressure of the dust fluid instead of any dark energy component. 
This fluid represents an unified description of the dark sector plus the baryon component 
in a similar way than the Chaplygin gas model (see for instance [291 PfO] and references 
therein). One of the advantages of the model is that it solves automatically the cosmic 
coincidence problem because there is not any dark energy component. The explicit 
form of the bulk viscosity has to be assumed a priori or obtained from a known or 
proposed physical mechanism. In the present work we choose the first option to explore 
the possibilities, but with the idea to explore the second option in future works. 

Different aspects of this model has been also studied in [2"?t HI] an d references 
therein. There is already a large amount of works on viscous fluids in cosmology 
with different parametrizations and motivations. For instance, the particular case of 
a constant bulk viscosity (the simplest parametrization) it has been carefully studied in 
detail recently [I2J [43J . Some other works with the constant parametrization has been 
also analyzed in [21], [H]-[17]. Other parametrizations or more general approaches for 
the viscous models have been also proposed and studied at [2~5]-[3T)]. [IE]- [5]]. Another 
general approach to the bulk viscous cosmologies called "fluids with inhomogeneous 
equation of state" have been also proposed and studied by [52] . 

In section [2] we present the main formalism of bulk viscous fluids in General 
Relativity (GR) and we apply it to the bulk viscous matter-dominated model studied 
in this work. We write the expression of H in terms of the density parameter of the 
viscous matter component in function of the redshift. With the expression of H we find 
the explicit expression for the scale factor that will be used in section [3] to analyze the 
possible scenarios for the Universe according to the values of the dimensionless viscous 
coefficients Co and Ci- 

In sections H] and [5] we analyze the curvature scalars and the matter density to get 
a better understanding of the behavior of the Universe according to the present model 
when it is evaluated at the best estimated values of Co and Ci- In section [6] we do a 
brief review of the local second law of thermodynamics and its possible violation for 
these models. Section [JJ presents the SNe la probe to be used to constrain the model 
and to compute the best estimated values of (Co? Ci) • Finally, in section [8] we give our 
conclusions. 
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2. Cosmological model of bulk viscous matter-dominated Universes. 

We analyze the cosmological model in a spatially flat Universe framework. For that, we 
use the Weinberg formalism [53] of imperfect fluids. So, in the present work we consider 
a bulk viscous fluid as source of matter in the Einstein fields equations G^ v = 8-nGT^, 
where G is the Newton gravitational constant. 

The energy-momentum tensor of the bulk viscous matter component is that of an 
imperfect fluid with a first-order deviation from the thermodynamic equilibrium. It can 
be expressed as [3U [53] : 

= Pm U^Uv + (g^ + UpU^Pn , (1) 

where 

P^ = P m - CV„^ , (2) 

m m is the four-velocity vector of an observer who measures the energy density p m and 
P m is the pressure of the fluid of matter, the term g^ v is the metric tensor, the subscript 
"m" stands for "matter" component. The term £ is a bulk viscous coefficient that arises 
in a fluid when it is out of the local thermodynamic equilibrium and that induces a 
viscous pressure equals to — CV u u u [311 E7]-[39]. The term P^ is an effective pressure 
composed by the pressure P m of fluid plus the bulk viscous pressure. 

The effective pressure (E]) was proposed by Eckart [54J in 1940 for relativistic 
dissipative processes in thermodynamics systems out of local equilibrium. Later, Landau 
Sz Lifshitz developed also an equivalent formulation [55J. 

The Eckart theory has some problems in its formulation, for example, all the 
equilibrium states in this theory are unstable [56J, other issue is that signals can 
propagate through the fluids faster than the speed of light [HTJ [58]. To correct the 
problems of the Eckart theory, Israel-Stewart [59|l60] developed in 1979 a more consistent 
and general theory that avoids these issues from which the Eckart theory is the first- 
order limit when the relaxation time goes to zero. In this limit the Eckart theory is a 
good approximation to the Israel- Stewart theory. 

In spite of the problems of the Eckart theory, but taking advantage of the 
equivalence of both theories at this limit, it has been widely used by several authors 
because it is simpler to work with this than with the Israel-Stewart one. In particular, 
it has been used to model bulk viscous dark fluids as responsible of the observed 
acceleration of the Universe assuming that the approximation of vanishing relaxation 
time is valid for this purpose (see, for instance [25]- [27], [29| |30j |32| |i2j |33] ) . 

It is important to point out that Hiscock et al [21] showed that a flat Friedmann- 
Robertson- Walker cosmological model with a bulk viscous Boltzmann gas expands faster 
when the Eckart framework is used than with the Israel-Stewart one. He suggests that 
the inflationary acceleration due to the bulk viscosity could be just an effect of using 
the uncausal Eckart theory and not a real acceleration. However, posterior studies seem 
to show that this conclusion could be wrong because consistent inflationary solutions 
have been found using the Israel- Stewart theory [23J. 
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The model could have a problem related with the cosmological density 
perturbations, where a rapid late time decay or irregular perturbations arise, leading 
to large modifications of the model Lambda Cold Dark Matter (ACDM) predictions 
for the matter power spectrum, CMB and weak lensing tests, as shown by B. Li et 
al. [49J. However, Hipolito et al. [51] have shown that it is possible to have models 
with bulk viscosity proportional to ( oc p u that are compatible with the data from the 
2dFGRS and the SDSS surveys. Due to the problem with the density perturbations, in 
the present work we assume that the viscosity appears only till recent times. 

It is convenient to mention that for irreversible processes, D. Pavon et al. [61] 
have developed a more general formulation than the Israel- Stewart theory, where the 
temperature could not be associated to the thermal equilibrium. However, in the 
cosmological scenario, this general formulation has not been totally explored. 

Let's start assuming a spatially flat geometry for the Friedmann- Robert son- Walker 
(FRW) cosmology as suggested by WMAP [62] 

ds 2 = -dt 2 + a 2 (t) (dr 2 + r 2 dVl 2 ) (3) 

where a(t) is the scale factor. Moreover, the conservation equation for the viscous fluid 
is 

u v V vPm + (p m + P* ) V v u v = (4) 

We write this conservation equation using the FRW metric as 

p m + 3P(p m + P* ) = (5) 

where H = a/a is the Hubble parameter. The dot over p m and a means time derivative. 
For the matter component, we set P m = 0, in concordance with a pressureless fluid. The 
bulk viscous pressure —CV u u u can be written as —3(H. So, the conservation equation 
fl5]) for the pressureless viscous fluid becomes 

p m + 3P(p m - 3CP) = (6) 

The Friedmann equations for this model are 

H = —P™ ( 7 ) 

d AixG , ^ >TT x ,„ s 
- = ^(p m -9CP) (8) 

Using the ansatz C = Co + Ci-^ an d the first Friedmann equation ([7]), we can write the 
term (p m — 3£P) of equation (jBJ) as 

p m - 3CP = (l - |j p m - (24^1/2 CoPm /2 (9) 
where we have defined the dimensionless bulk viscous coefficients Co an d Ci as 

Co=(^^)Co, Ci = (24ttG)Ci (10) 
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and H is the Hubble constant. Using the expression we can rewrite the conservation 
equation (jfJJ) in terms of the scale factor as 



dpm . (3 - Ci) 
da a 



3-^p Co 
(24ttG') 1 /2 a 



P 



*/ 2 = 



'11' 



We divide the ordinary differential equation (ODE) ( fill by the critical density today 
p° crit = 3H$/8nG yielding 



da 



+ (3-Ci)fi m -Co^ /2 = 



(12) 



where we have defined the dimensionless density parameter Q m = p m / P^if Using the 
relation between the scale factor a and the redshift l + z = l/awe can write the ODE 
AHD as 



;i + ^)^ + (Ci-3)fi m + Co^ /2 = 



(13) 



where the function to be calculated in this ODE is Cl m (z). This ODE has the exact 
solution 

2 



Co 



3-Ci 



(l + s) 



(3-C0/2 



+ 



Co 



3-d 



for Ci ^ 3 (14) 



where we have used the initial condition f2 m0 ee Q m (z = 0) = 1 |jj in concordance with a 
flat Universe. 

On the other hand, we divide the first Friedmann equation by the critical density 
today and then substitute there the expression ffl4|) to obtain 



3-Co-Ci (l + ^) (3 - Cl)/2 + Co 
H(z) = H ± ; 



3-Ci 



Or in terms of the scale factor 



H(a) = H n 



1 - 



3-Ci 



a (Ci-3)/2 + Co 



3-Ci 



for Ci ^ 3 



for Ci ^ 3 



(15) 



(16) 



% In the present work we assume tt m o — 1 like many other authors (see for instance [63l [64] and 
references therein). To be consistent with the dynamical observations measured in clusters of galaxies 
that suggest il m o ~ 0.26, it is assumed that this amount of matter is bounded in structures like galaxies 
and cluster of galaxies and that the remaining matter f2 m o ~ 0.74 is distributed homogeneously in the 
flat space. 
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3. The scale factor. 



We analyze the predicted behavior of the scale factor and the evolution of the Universe 
according to the present model. So, from the definition of the Hubble parameter 
H(a) = a/a we write the expression <HM as 



Ida 

a at 



1 



Co 



a 



(Ci-3)/2 



Co 



3 - Ci J 3 - Ci 

We define 7 = Co/ (3 — Ci) an d P = Ci — 3, and integrate the equation ffTTl) as follows 

da' rt 



(17) 



f H dt f = H (t - t c 



I, a' [(1-7)0^/2 + 7] Jt0 
where t is the cosmic time. Solving this integral we obtain 

7 



11 



7/3 



7 + 



# (* - h) 



(19) 



Solving for a(t) and substituting the values of 7 and we arrive to 

2/(3-00 



a(t) 



Co + Ci - 3 + (3 - Ci) exp k H (t - t ) 



Co 



; for Co t^O, Ci^3 (20) 



Also, we compute the first and second derivatives of expression ( 120]) with respect 
to x = H (t — t ) to study the accelerated or decelerated epochs of the Universe and 
the transitions between them 



da 
dx 



Co + Ci - 3 + (3 - Ci)e' 



Cox/2 



Co 



(Ci-l)/(3-Ci) 



d 2 a 
dx 2 



-e 
2 



0W2 



Co + Ci - 3 + 2e^ 2 



Co + Ci - 3 + (3 - C i)e 
Co 



Cox/2 



-, 2(Ci-2)/(3-C0 



(21) 



(22) 



We have analyzed the behavior of the scale factor in all the possible combinations 
of values for (Co, CO- We have found that from all the possibilities, the cases where there 
is a Big-Bang as the origin of the Universe, followed by a decelerated expansion epoch at 
early times and with a transition to an accelerated epoch in recent times, corresponds to 
the cases with values (Co > 0, Ci < 0) an d (0 < Co < 3, Ci = 0). We have realized also 
that the best estimated values for (Co, Ci), or Co alone (when we set £1 = 0) computed 
using the latest SNe la data set, lie precisely in these two ranges of values respectively 
[see table [2]. We present below the details of the analysis of these two cases. 

It is important to highlight that in general for any other case different of the two 
cases (Co > 0,Ci < 0) (with Co + Ci < 3), and (0 < Co < 3,Ci = 0), the Universe does 
not have the same behavior, i.e., for some other cases there is not a Big-Bang or a late 
time accelerated epoch with a early time deceleration period, etc. 




Figure 1. Plot of the scale factor with respect to Ho(t — to) for different values of 
(CojCi) when ( > and (1 < but with the condition Co + Ci < 3 [see expression 
(1201) ]. The model predicts an eternally expanding Universe that begins with a Big- 
Bang followed by a decreasing decelerated expansion until a time when the deceleration 
vanishes and then there is a transition to an accelerated expansion epoch that is going 
to continue forever, so that a — > oo when t — > oo. See section [3] for details. The 
points indicate the time when the transition between the decelerated expansion epoch 
to an accelerated one happens [see expression ()24j) ]. The dashed line corresponds to a 
de Sitter Universe that is obtained when Co + Ci = 3. 



First, we describe the case (Co > 0, Ci < 0). There are two kinds of behaviors 
depending on the inequality: 

(a) Co + Ci < 3, (b) Co + Ci > 3 (23) 

Figures HH5] show the behavior of the scale factor for these subcases using different 
values for (Co, Ci)- 

For the subcase (a) the model predicts an eternally expanding Universe beginning 
with a Big-Bang^ in the past followed by a decelerated expansion where the deceleration 
is decreasing until a time when it vanishes and then there is a smooth transition to 
an accelerated expansion epoch that is going to continue forever, so that a — > oo when 
t — t- oo (see figures dJUJ) . 

The cosmic time of transition "ttrans" between the decelerated to the accelerated 
expansion epochs can be calculated equating to zero the expression ff22l yielding 

t tr ans = t + -^luf 3 "^~ C ^ (24) 

-noCo \ / 

§ In the present work we assume as a Big-Bang the state of the Universe where the scale factor is zero 
in some cosmic time teB in the past of the Universe, i.e., a(ijiB ) = 0, where in addition, the curvature 
scalar R and the matter density p are singulars, i.e., R(a = 0) = p(a = 0) = oo (see sections [4] and [5] 
for details). The subscript "BB" stands for "Big-Bang". 
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Figure 2. Plot of the first and second derivative of the scale factor with respect to 
x = Ho(t — to) for different values of (Co,Ci) when £ > and Ci < but with the 
condition Co + Ci < 3 [see equations ([2T]) and ([22])]. These plots show the decelerated 
and accelerated expanding epochs as well as the transition between them (see section 
[3]). For the left figure, the minimum value of each plot corresponds to the cosmic 
time when the transition happens, that for the right figure corresponds when a plot 
intercepts the x-axis. 



To find the expression of the transition between the decelerated to the accelerated 
expansion epochs in terms of the scale factor we derive with respect to a the expression 
([IE]) and using the fact of that H = a/a we obtain 



da 
da 



Ci-i 



1 



Co 



3-Ci 



a (Ci-3)/2 + 



Co 



3-Ci 



(25) 



We equal to zero the equation ([25]) yielding 

2/(Ci-3) 



a t 



2Co 



(Ci-l)(Co + Ci-3) 



(26) 



The transition happens in the past of the Universe if 1 < Co + Ci < 3, in the future if 
Co + Ci < 1 an d today if Co + Ci = 1- 

On the other hand, using a(t^) = 0, we find the cosmic time £bb when the Big- 
Bang happens 



t 



BB 



t + 



H Q 



^ In 1 - 



Co 



(27) 



o 



3-Ci 

We define the Age of the Universe as the elapsed time between the time £bb until 
the present time t . So, using the expression (1271) . the age of the Universe is given by 



Age = \t - t 



BB | 



■^ln 1 



Co 



#oCo 



3-Ci 



(2f 



The best estimated values for the bulk viscous dimensionless coefficients (Co, Ci) 
computed using the type la SNe data correspond to this subcase (see tableland section 
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Figure 3. Plot of the scale factor a(t) with respect to Hg(t — to) [see expression (|20l) ] 
when it is evaluated at the best estimated values of Co and Ci (see table [2]). The black 
line corresponds to the values (Co = 8.58 ± 1.2, Ci — —5.96 ± 1.1) coming from the 
use of a Dirac Delta prior distribution to marginalize over the Hubble constant H$. 
The red line corresponds to the values (Co = 4.38 ±1.5, £i — — 2.16 ± 1.4) coming 
from the use of a Constant prior distribution to marginalize over Hq (see table [2]) . 
The plot shows how the model, evaluated at the best estimates, predicts a Big-Bang 
as origin of the Universe, followed by a decelerated expanding epoch at early times, 
with a transition to an accelerated epoch at late times that will continue forever. The 
points indicate the cosmic time when the transition between the decelerated epoch to 
an accelerated one happens [see expression (|24|) and table Q]. The grey and pink bands 
correspond to evaluate the expression (|20|) at the errors in the estimations for bulk 
viscous coefficients, i.e., for the Dirac delta case (pink band) the expression (|20|) was 
evaluated at (Co = 8.58 + 1.2, Ci = —5.96 + 1.1) to plot the upper limit of the band 
and (Co = 8.58 — 1.2, Ci = —5.96 — 1.1) for the lower one. In the same way, for the 
constant prior (grey band) (Co = 4.38 ± 1.5, Ci = —2.16 ± 1.4) for the upper and lower 
limits respectively. 



[7] for details). The figures [3] and H] show how the model, evaluated at the best estimates, 
predicts a Big-Bang as origin of the Universe, followed by a decelerated expanding epoch 
at early times, with a transition to an accelerated epoch at late times that will continue 
forever. 

Tabled] shows the values for the age of the Universe, evaluated at the best estimates, 
as well as the value of the redshift of transition between deceleration-acceleration epochs. 
From table [1] we can see that for the case of both (Co, Ci) as f ree parameters (i.e., when 
Ci 7^ 0), and for both marginalizations over the Hubble constant H , the estimated 
values for the age of the Universe (11.72 and 10.03 Gyr) are smaller than the mean 
value of 12.9 Gigayears coming from the oldest globular cluster observations but still 
inside of the confidence interval of these: 12.9 ±2.9 Gigayears (see section [7] and table [1] 
for details). 

For the subcase (b) the model predicts an eternally expanding Universe and this 
expansion is always accelerated along the whole history of the Universe (see figure |5]). 
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Figure 4. Plot of the first and second derivative of the scale factor a(t) with respect 
to x = Ho(t — to) [see expressions (|2Tj) and (|22|) ] when they are evaluated at the 
best estimated values of Co and Ci (see table [2]). The points indicate the time when 
the transition between the decelerated expanding epoch to an accelerated one happens 
[see expression (l24l) and table Q]. For the left figure, the minimum value of each 
plot corresponds to the cosmic time when the transition happens, that for the right 
figure corresponds when a plot intercepts the x-axis. The legends " 'delta prior", 
11 constant prior" and 11 constant prior, Ci — 0" label the different marginalizations 
used to compute the best estimated values (see table [2]). 



In this subcase there is not a Big-Bang as the origin of the Universe, instead of that, 
the Universe has a minimum value of the scale factor given by 



where the subscript "min" stands for "minimum" . From this minimum value a m i n the 
scale factor increases along the time until a — > oo when t — > oo. For this subcase the 
age of the Universe is not defined [see equation (125)) ]. 

Now, for the second case (0 < Co < 3, Ci = 0)> ft nas a l so the same behavior 
than the case Co > 0, Ci < (with Co + Ci < 3), it is, the model predicts a Big-Bang 
as origin of the Universe, followed by an early time decelerated expanding epoch, with 
a transition to an accelerating epoch that is going to continue forever (see figures ID 
and |7]) . When we set Ci = 0, the best computed value for the other viscous coefficient 
is Co = 1-9836 ± 0.066 (see table [2j. The transition between deceleration-acceleration 
happens when the scale factor has a value of a t = 0.4 [see table [Hand expression (1261)]. 
In this case best estimated value for Co predicts an age for the Universe of 14.82 Gyr, 
that is in perfect agreement with the constraints of globular clusters. An important 
feature of this case is that it does not violate the second law of the thermodynamics 
(see section [6]). This particular case has been also studied in [42] . 




(29) 
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Figure 5. Plot of the scale factor with respect to Ho(t — to) for different values 
°f (CoiCl) when ( > and £i < but with the condition of (o + Ci > 3 [see 
equation (1201) ]. The model predicts an eternally expanding Universe and this expansion 
is always accelerated. There is not a Big-Bang as the origin of the Universe, 
instead of that, the Universe has a minimum value of the scale factor given by 
Omin = [(Co + Ci — 3) / Co} 2 ^ 3 ^^ 1 ^ when t — > — oo. From this minimum value a m ; n 
the scale factor increases along the time until a — > oo when t — > oo. See section [3] for 
details. 



Viscous model £ = Co + Ci-^ 



Prior over Hq 


Age Universe 


a t z t 


Constant 
Dirac delta 
Constant (Ci = 0) 


11.72 Gyr 
10.03 Gyr 
14.82 Gyr 


0.61 0.63 
0.65 0.52 
0.40 1.47 



Table 1. Summary of the computed central values of the age of the Universe, the scale 
factor of transition between deceleration-acceleration epochs "at" and its corresponding 
value in redshift "z t " , coming from the best estimates of £o and £i [see table [2]. The 
first column shows the prior distribution assumed for Hq to marginalize. The second 
column shows the estimated central values for the age of the Universe [see expression 
(|28|) ] given in units of Gigayears (Gyr). These magnitudes were obtained assuming a 
value of Hq = 72 (km/s)Mpc~ 1 in the expression (|28|) as suggested by the observations 
of the Hubble Space Telescope (HST) [70], assuming a value of a year of 31556925.2 
seconds (a tropical year) and a megaparsec = 3.0856776 x 10 19 km. The third column 
shows the estimated value for scale factor at when the transition between the decelerated 
expanding epoch to the accelerated expanding one happens [see expression (|26p]. The 
fourth column shows also the value of this transition expressed in terms of the redshift 
and labeled as z t (the relation is given by 1 + z = 1/a). 
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Age of the universe 




100 



Figure 6. Plot of the age of the Universe in terms of the Hubble constant H when 
it is evaluated at the best estimated values for (CchCi) [ see expression (|2"8|) and tables 
[T]and |2j. Hq is given in units of (km/s)Mpc _1 and the age in units of Gigayears 
(Gyr). It was assumed a value of a year of 31556925.2 seconds (a tropical year) and 
a megaparsec = 3.0856776 x 10 19 km. The points indicate the computed value for 
the age of the Universe when it is assumed a value of Hq = 72 (km/s)Mpc~ 1 , as 
suggested by the observations of the Hubble Space Telescope (HST) [70]. The vertical 
lines correspond to the interval H = [55, 75] (km/s)Mpc , it is the permitted region 
according to values of H consistent with the distance moduli used to derive ages for 
Galactic globular clusters from the Hipparcos parallaxes (see ref. [66]). The horizontal 
lines correspond to the constraint for the age of the Universe from the oldest globular 
clusters (Age= 12.9 ± 2.9 Gyr [66 ). So, the shaded area is the consistent region for 
the age of the Universe. Table [T] summarizes the estimated values for the age of the 
Universe. 



4. The curvature 

To study the possible singularities of the model we analyze the behavior of the curvature 
scalar R, as well as the contractions of the Riemann tensor R a ^ s R a ^ and the Ricci 
tensor R al3 R a p. We find that there is not any particular singularity other than that 
of the Big-Bang, that could indicate possible issues for the model coming from the 
curvature. About Weyl tensor, it turns out that C a ^ s C a p 7 s = in general for the FRW 
metric. The details of the analysis is shown below. 

4-1. The curvature scalar R 

The curvature scalar R = R a a for a spatially flat Universe is given as 

a 



R = 6 



- + H 



(30) 
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Figure 7. Plot of the scale factor a(t) with respect to Hg(t — to) [see expression (|20l) ] 
evaluated at the best estimated value Co = 1.98 ± 0.06 (see table[2|). It corresponds to 
the case when £i is set to zero and it is computed the best estimate for £ alone, using 
a constant prior to marginalize over Hq. The legend "Constant prior, £i = 0" indicates 
this. The plot shows how the scale factor, evaluated at the best estimated value for Co 
(with £i =0), predicts a decelerated expanding epoch just after de Big-Bang followed 
by an accelerating epoch that is going to continue forever. The transition between 
deceleration-acceleration happens when the scale factor has a value of a t — 0.4 [see 
table [Hand expression (|26|) ]. it is indicated with the black point in the plot. The thin 
pink band corresponds to evaluate the expression (|20|) at the statistical error of the 
estimation of Co- 



The term d/a corresponds to the second Friedmann equation, i.e., 
a A-kG 



(31) 



a 3 

We substitute the expressions for: P* a = —3(H [see eq. ([2}], C = Co + Ci^ an d 
Pm — {3/8ttG)H 2 from the first Friedmann equation at the expression (l3"Tj) to obtain 



a 



;i2vrG)(| + Ci 



(32) 



Using the expressions of the dimensionless bulk viscous coefficients Co an d Ci [ see 
expressions (flOT) ] yield 

" ' (33) 



Co + Ci 



H 



We substitute the equation ( 1331) at ( 1301) to obtain 



R = 3H 



#oCo + #(Cl + 1; 



(34) 



Now, using the equation (Tl6l) for the Hubble parameter we arrive to 



R(a) 



3# 2 



(3-Ci) 2 
for Ci + 3. 



4Co + (Co + Ci - 3)a 



(Ci-3)/2 



(Ci + l)(Co + Ci - 3)a^" 3 )/ 2 - Co(Ci + 5)1 ) (35) 
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For the case C > 0, Ci < 0, when a ->• oo then R -> 12[if C /(3 - Ci)] 2 - When 
in addition Co + Ci < 3, then the curvature scalar diverges when the scale factor a is 
equal to zero. This supports the existence of a Big-Bang in the past of the Universe 
(see figure [8] and section . 

However, when Co + Ci > 3 (with Co > 0, Ci < 0) there is not Big-Bang in the past 
of the Universe because the scale factor never becomes zero in the past, the minimum 
value that the scale factor can reach is a min = [1 + (Ci — 3)/Co) ^ when t — >■ — oo. 
So, the curvature scalar does not diverge but it is zero when the scale factor is equal to 
a min . Then, there is not the conditions for a Big-Bang when Co + Ci > 3 (see figured]). 

When Co + Ci = 3 the curvature scalar has the constant value of V2Hq along the 
whole evolution of the Universe. When Ci > —1 the curvature scalar is positive as a — > 
and negative when Ci < — 1- When Ci = ~ 1 then the curvature scalar is positive as 
a — > if Co < 4 and negative if Co > 4. The present central value of the curvature scalar 
is R/Hq = 9.66 and R/Hq = 10.87 when constant and Dirac delta priors for H are 
used, respectively. 

For the other important case (0 < Co < 3, Ci =0), we find that the curvature scalar 
diverges when a — > (the Big-Bang singularity). From the infinite value of R(a = 0) it 
decreases as the Universe expands until to reach its minimum value of R = f-^oCo when 
a — > oo (see figured]). 

4.2. The scalar R a ^ & R aM& 

We analyze the behavior of the contraction of the Riemann tensor with itself, i.e., 
K = R a/3lS R a p^s ■ In general, for a FRW metric in a curved space-time K is found to be 



where k characterizes the constant spatial curvature of the FRW metric; for our model 
k = 0. So, using the fact of that H = (a/ a) and that a /a corresponds to the second 
Friedmann equation given by expression ( l33l) we rewrite the equation ( l36l) as 




(36) 




(37) 



where H is given by the expression (fl6l) . that we substitute to obtain 




(38) 



For the case (Co > 0, Ci < 0), we find that K diverges to infinity when a — > and 
when a — > oo then K — > 24[_ff Co/(Ci — 3)] 4 . For the best estimates, the present central 
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2.0 

Figure 8. Plot of the curvature scalar R with respect to the scale factor for different 
values of (Co, Ci) f° r the cases (Co > 0, Ci < 0) and (0 < Co < 3, Ci = 0) [see equation 
([55]) ] . The short and long dashed lines correspond to evaluate R at the best estimates 
for (Co, Ci) when the constant and Dirac delta priors are used to marginalize over Hq, 
respectively. The solid black line corresponds to evaluate R at the best estimate of 
Co alone, when we set Ci = and marginalize Hq using a constant prior. In general, 
the curvature scalar diverges when the scale factor is equal to zero for the subcase 
Co + Ci < 3 (with Co > 0, Ci < 0). For the subcase Co + Ci > 3 (with Co > 0, Ci < 0) the 
minimum value that the scale factor reaches is a m i n = [1 + (Ci — 3)/Co] 2 ^ 3_< ° 1 - ) when 
t — > — oo. So, for this subcase the curvature scalar does not diverge but it is zero when 
the scale factor is equal to a m j n . Then, there is not the conditions for a Big-Bang. 
When Co + Ci =3 the curvature scalar has the constant value of 12Hq along the whole 
evolution of the Universe. When Ci > — 1 the curvature scalar is positive as a — » 
and negative when Ci < — 1- When Ci = — 1 then the curvature scalar is positive as 
a -> if Co < 4 and negative if Co > 4. When a ->■ oo then R ->■ 12[iJ Co/(3 - Ci)] 2 - 
Finally, for the case (0 < Co < 3, Ci = 0), the curvature scalar diverges when a = 
(the Big-Bang singularity), and when a — > oo then R = l-ff^Co- 

value for the Riemann contraction is K/Hq = 16.48 and K/Hq = 19.91 when constant 
and Dirac delta priors are used, respectively. Figure M shows the behavior of K with 
respect to the scale factor, evaluated at the best estimates. 

For the other case of interest (0 < Co < 3, Ci = 0), we find that the scalar K 
is infinite when a — > (the Big-Bang singularity), it decreases from K(a = 0) as the 
Universe expands until to reach its minimum value of K — (8/27)(ffoCo) 4 when a — > oo 
(see figure [9]). 

4.3. The scalar R a/3 R a p 

The scalar derived of contracting the Ricci tensor with itself, A = R a/3 R a /3, for a curved 
FRW metric has the form 

A = R afS R a 8 = 12 
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+ —(a 2 + k) + 



+ 2k 



+ 



(39) 
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Figure 9. Plot in logarithmic scale of the scalars of curvature K (left) and A (right), 
coming from the contraction of the Riemann tensor K = R a P~< & R a j3 7 s, and the Ricci 
tensor A = R a/3 R a p, with respect to the scale factor [see expressions and (ptTj) 
respectively]. Both scalars are evaluated at the best estimates for (CojCi) together, 
and Co alone (when we set Ci =0)- See tabled Note that in both figures and for the 
six plots, K and A diverge when a — > (the Big-Bang singularity) and when a — > oo 
then K and A go to a constant minimum value. 



For our model k = 0. We substitute a/a by H and a/ a by the expression ( )33|) to obtain 



A = 3H 2 



+ 2tf CoCi# + (3 + Cf)H 2 (40) 



We put the expression for the Hubble constant ( fTBl in f j40~l) yielding 
A(a) = 3H i A - {Co« 3/2 + (3 - Co - Ci)a &/a V x (41) 



(Ci - 3) 4 a 6 

x {l2C 2 a 3 + 6Co(l + &)(3 - Co - C> (3 ^ l)/2 - (3 + Ci 2 )(3 - Co - Ci) 2 ^ 1 } 

For the case (Co > 0, Ci < 0), when a — > then A — > oo. And when a — > oo then 
A — > 36[-£/oCo/(Ci — 3)] 4 . For the best estimates, the present central value of the Ricci 
contraction is A/Hq = 23.81 and A/Hq = 29.66 when constant and Dirac delta priors 
are used, respectively. Figure M shows the behavior of A with respect to the scale factor. 

For the other case (0 < Co < 3, Ci = 0), the scalar A is infinite when a — > (the 
Big-Bang singularity), it decreases from A(a = 0) as the Universe expands until to reach 
its minimum value of A = (4/9)(H ( ) 4 when a — > oo (see figure |9]). 



5. Matter density p m 

We analyze also the equation of the matter density to have a better understanding of the 
behavior of the Universe according to the model. From the first Friedmann equation ([7]) 
we have 

Pm(«) = (g^) H\a) (42) 
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Figure 10. Plot of the matter density p m with respect to the scale factor for different 
values of (Co, Ci) where Co > and Ci < (left panel) and also in the best estimates for 
(Coi CO (right panel) [see equation (|43| and table [2]. For the subcase Co + Ci < 3 (with 
Co > 0, Ci < 0) the matter density diverges when the scale factor is equal to zero. This 
supports the assumption of a Big-Bang as the beginning of the Universe (see section 
[3]). However, for the subcase Co + Ci > 3 the minimum value that the scale factor 
reaches is a min = [(Co + Ci — 3)/Co] 2 ^ 3- ^ when t — > — oo. The matter density is zero 
when the scale factor is equal to a m i n - Then, there is not the conditions for a Big-Bang 
in this subcase. In both subcases, when a — > oo then p m = /9° rit [Co/(3 — Ci)] 2 - F° r the 
subcase Co + Ci = 3, the matter density has the constant value of Pcrif F° r the right 
panel, the short and long dashed lines correspond to evaluate the matter density at the 
best estimates when it is used the constant and Dirac delta priors to marginalize over 
Hq, respectively. The solid black line corresponds to the best estimate for Co> setting 
Ci = and using a constant prior to marginalize over Hq. 

Using expression for the Hubble parameter ( IT6|) at the expression above we arrive to 

2 

; for Ci ^3 (43) 

When Co > 0, Ci < 0: For the subcase Co + Ci < 3 (with Co > 0, Ci < 0) the matter 
density diverges when the scale factor is equal to zero. This supports the assumption 
of a Big-Bang as the beginning of the Universe (see section [3] and figure fTOl) . 

However, for the subcase Co + Ci > 3 the minimum value that the scale factor reaches 
is a m in = [(Co + Ci — 3)/Co] 2 ^ 3 ^^^ when t — > — oo. The matter density is zero when the 
scale factor is equal to a m i n - Then, there is not the conditions for a Big-Bang in this 
subcase. In both subcases, when a — > oo then p m = Pc rit [Co/(3 — Ci)] 2 - 

For the subcase Co + Ci = 3, the matter density has the constant value of p° rit . 
Figure [10] shows the behavior of the matter density with respect to the scale factor. 

For the other case of interest, (0 < Co < 3, Ci = 0), the matter density is infinite 
when a — > (the Big-Bang singularity), from this infinite value p m decreases as the 
Universe expands until to reach its minimum value of p m = (Hq/24ttG)(q when a — > oo 
(see figure [10]) . 
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Figure 11. Plot of the total dimensionless bulk viscosity C = Co + (iH with respect to 
the scale factor (left panel) and the redshift (right panel) [see equation (|48|)] when it is 
evaluated at the best estimated values for (Co,Ci) ( see table [2]). At the early Universe 
(small values of a, or greater values of z) the total bulk viscosity C is negative for the 
case of "delta prior" and "constant prior" . Negative values for the total bulk viscosity 
C implies a violation of the local second law of thermodynamics at early times [see 
equation (l47l) ]. So, the model, ( = Co + CiH, works well only for recent times of the 
Universe. The plot "constant prior, Ci = 0" corresponds to the model C — Co evaluated 
at the best estimated value for Co (we have set Ci =0)- This case does not violate 
the local second law of thermodynamics, the supernova data indicates a positive value 
for Co and it is constant along the whole history of the Universe. The transition to 
a positive values in recent times of the model C = Co + Ci-^ happens when the value 
of the scale factor a np is given by equation (j49|) . And when a — > oo, the total bulk 
viscosity for this case is C = 3Co/(3 — Ci)- 



6. Thermodynamics and the local entropy. 

The law of generation of local entropy in a fluid on a FRW space-time can be written 
as [53] 

T V u s u = C(V v u u ) 2 = 9H 2 ( (44) 

where T is the temperature and V u s u is the rate of entropy production in a unit volume. 
With this, the second law of the thermodynamics can be written as 

TV u s u > (45) 
so, from the expression (jHJ), it implies that 

C > (46) 
Thus, for the present model the inequality ()46l) can be written as 

C = Co + CiH > (47) 

Using the expression for the Hubble parameter ([TBI) we find that the expression for the 
total bulk viscosity ((a) is given as 
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where we have defined the total dimensionless bulk viscous coefficient ( = (24nG/ H )(. 
Figure [11] shows the behavior of C with respect to the scale factor when (Co? Ci) are 
evaluated at the best estimated value (see section [7] and table [2]). We find that for 
the best estimates of (Co, CO? t ne total bulk viscosity is negative in early times of the 
Universe and positive at late times. For the best estimate of Co alone (when we set 
Ci =0), the total bulk viscosity is positive and constant along the cosmic time, because 
in this case C = Co- 

According to equation (I4"T1) . negative values for the total bulk viscosity C implies 
a violation of the local second law of thermodynamics. So, the model C — Co + Ci-^, 
evaluated at the best estimates, violates the entropy law at early times, but the model 
C = Co? evaluated at the best estimate for Co (with Ci = 0), does not. 

The transition between negative to positive values of the total bulk viscosity for 
the case C = Co + Ci-^ happens when the scale factor has the value 

2/(6-3) 

(49) 

The subscript "np" stands for "negative to positive" values. For the best estimates 
for (Co, Ci) using the SNe la probe we find that the value of the scale factor a np when 
transition happens are a np = 0.45 and a np = 0.57 when the constant and Dirac delta 
priors are assumed to marginalize over Ho, respectively. When a — > oo then the total 
bulk viscosity is C = 3Co/ (3 — Ci)- 



'-up 



3Co 



Ci(Co + Ci-3) 



7. Type la Supernovae test. 

We test and constrain the viability of the model using the type la Supernovae (SNe 
la) observations. So, we calculate the best estimated values for the parameters Co 
and Ci an d the goodness- of- fit of the model to the data by a ^-minimization and 
then compute the confidence intervals for (Co, Ci) to constrain their possible values with 
levels of statistical confidence. We use the "Union2" SNe la data set (2010) from "The 
Supernova Cosmology Project" (SCP) composed by 557 type la supernovae [65] . 

We use the definition of luminosity distance d^ (see [3j [5], [67]-[69]) in a flat 
cosmology, 

~ ~ r z dz' 

d L (z, Co, Ci, H ) = c(l + z) / u( , ~ ~, (50) 

Jo H(z', Co, Ci> Ho) 

where H(z, Co, Ci, -^o) is the Hubble parameter, i.e., the expression f[T5"j) and 'c' is the 
speed of light given in units of km/sec. The theoretical distance moduli for the k-th 
supernova with redshift z k is defined as 

f^(z k , Co, Cx, H o) = rn - M = 5 log 10 



djXzk, Co, Ci, Hp 
Mpc 



+ 25 



(51) 
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Figure 12. Confidence intervals for the dimensionless coefficient (CojCi) °f a bulk 
viscous matter-dominated Universe model with bulk viscosity of the form C = Co + Ci-^- 
We used the SCP (2010) "Union2" SNe data set composed of 557 type la SNe [65], 
to compute the best estimates Co = 4.389 ± 1.56 and Ci — —2.166 ± 1.42; they are 
indicated with the blue point in the figure. The confidence intervals shown correspond 
to 68.3%, 95.4% and 99.73% of probability. The grey region (upper region) corresponds 
to the values for (Co,Ci) to which the model ( = (o + (iH predicts that there is not 
a Big-Bang in the past, it is, the scale factor never becomes zero in the past of the 
Universe (see figure [5]) , as well as the curvature scalar R did not diverge in the past 
(see figure [5]). The boundary line between the grey and white regions corresponds to a 
de Sitter Universe, it is obtained when Co + Ci = 3. The white region (central region) 
indicates the values to which the Universe is accelerating today, and the green region 
(lower) corresponds to deceleration today. Note that the confidence intervals of the 
best estimates for (Co, Ci) ne precisely in the accelerating region today with 99.73% of 
probability, at least. The dashed blue lines show two values for the age of the Universe 
(10 and 20 Gyr). The computed central value for the age using the best estimates for 
(Co, Ci) is 11.72 Gyr (see table [J). The Hubble constant H was marginalized assuming 
a constant prior distribution. 
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Viscous model ( = ( + (iH 



Prior over Hq 


Co Ci 


Xmin Xd.o.f. 


Constant 
Dirac delta 
Constant ((j = 0) 


4.38929 ± 1.568 -2.16673 ± 1.421 
8.58935 ± 1.217 -5.96505 ± 1.159 
1.9835 ±0.066 


542.38 0.977 
561.15 1.011 
544.587 0.979 



Table 2. Summary of the best estimated values of the dimensionless coefficients £o 
and Ci for a bulk viscous matter-dominated model with a bulk viscosity of the form 
£ = £o+Ci-ff • The best estimates were computed minimizing a x 2 function of a Bayesian 
statistical analysis using the SCP "Union2" 2010 compilation data set composed by 
557 type la SNe |65j . The first column shows the prior distribution function assumed to 
marginalize over the nuisance parameter Hq and then compute the best estimates for £o 
and (i once Ho has been marginalized. There are several ways to marginalize over Ho, 
in the present work we choose two: (1) "Constant prior" means that we do not assume 
any particular value for Hq, therefore, its probability distribution function is constant. 
(2) "Dirac delta prior" means that we do assume a specific value for Ho (in particular 
we used Ho — 72 (km/s)Mpc _1 ), therefore, its probability distribution function has 
the form of a Dirac delta centered at Ho = 72 (km/s)Mpc~ 1 . In the third row, the 
label "Constant (£i = 0)" means that we set £i = and use a constant prior for Hq, 
to compute the best estimated value for £ alone. The second column corresponds to 
the best estimated values computed for (Co,Ci), or Co alone (third row). The third 
column shows Xmin, the minimum value obtained for the x 2 function (|52|) . and also the 
X 2 function by degrees of freedom: Xdoi ■ ^ is defined as Xd o f = Xmin/ ( n ~P) where 
n is the number of data (in our case n — 557) and p is the number of free parameters 
of the model (for the first two raws p = 2 and for the third one p = 1). The errors are 
at 68.3% confidence level. Figures fT2l and fl3l show the confidence intervals. 



where m and M are the apparent and absolute magnitudes of the SNe la respectively, 
and the superscript 't' stands for "theoretical". We construct the statistical x 2 function 

as 

r ~ ~ i 2 

„ ~ - ™ \H fa, Co, Ci, Ho) - iik 

X 2 (Co, Ci, H ) = ~ 2 ~ (52) 

fc=i ° k 

where /if. is the observational distance moduli for the k-th supernova, o~\ is the variance 
of the measurement and n is the amount of supernova in the data set, in this case 
n = 557. 

Once constructed the x 2 function ( 1521) we numerically minimize it to compute the 
"best estimates" for the free parameters (Co, Ci) °f the model. The minimum value of the 
X 2 function gives the best estimated values of (Co, Ci) an d measures the goodness-of-fit 
of the model to data. The results are shown in table [2] and the confidence intervals for 
Co vs Ci are shown in figures [12] and [13j 

About the value of the Hubble constant Hq we marginalize it using two methods: 

Constant prior : We assume that Hq does not have any preferred value a priori, i.e, 
it has a constant prior probability distribution function. So, we solve analytically 
certain integrals using this assumption to obtain a x 2 -function like the expression 




Figure 13. Confidence intervals for the dimensionless coefficient (CojCi) °f a bulk 
viscous matter-dominated Universe model with bulk viscosity of the form £ = (q + ^iH. 
It is like the figure IT21 but now a Dirac delta prior was used to marginalize over Hq. 
The best estimates are £o = 8.58 ± 1.2 and £i = —5.96 ± 1.15 (see table [2]). Note how 
the confidence intervals of the best estimates lie precisely in the accelerating region 
today with 99.73% of probability. The computed central value for the age using the 
best estimates is 10.03 Gyr for this case (see table [T]). 



( 1521) but that it is not going to depend on the variable H anymore. This method 
is carefully described in the appendix A of [32"] . 

Dirac delta prior : We assume that H Q has a specific value (suggested by some 
other independent observation), so, its probability distribution function has the 
form of a Dirac delta, centered at the specific value. In particular we choose 
H = 72 (km/s)Mpc _1 as suggested by the observations of the Hubble Space 
Telescope (HST) [70] . In practice, it simply means to set H = 72 at the expression 
(|52|) [i.e., at expression (Tl5]) ] and with this H is not a free parameter of the model 
anymore. 



8. Conclusions. 

We have explored a bulk viscous matter-dominated cosmological model composed by a 
pressureless fluid (dust) with bulk viscosity of the form C = Co + CiH, representing both 
baryon and dark matter components. 




Figure 14. Plot of the normalized Probability Distribution Function (PDF) of Co for 
the case when we set Ci — 0. The computed best estimated value (the central value 
of the PDF) is Co = 1-9835 ± 0.066, the error is given at 68.3% confidence level (see 
tabled]). A constant prior was used to marginalize over H . 



One of the advantage of this model is that it resolves automatically the "Coincidence 
problem" because it is not introduced any dark energy component, i.e., it is enough to 
think of the Universe as composed by baryon and dark matter components as dominant 
ingredients of the Cosmos. 

The expansion of the Universe can be seen as a collection of states out of the thermal 
equilibrium during very short periods of time where it is produced local entropy. 

From all the possible scenarios predicted by the model according to different values 
of the dimensionless bulk viscous coefficients Co and Ci, we find two where the Universe 
begins with a Big-Bang, followed by a decelerated expansion epoch at early times, and 
with a transition to an accelerated expansion at recent times that is going to continue 
forever. These two scenarios are characterized for the ranges (Co > 0,Ci < 0) with 
Co + Ci < 3, and (0 < Co < 3, Ci = 0). This behavior of the Universe corresponds to 
what is actually observed and is close to the ACDM paradigm. 

On the other hand, we compute also the best estimated values for (Co, Ci) using the 
latest type la SNe data release. It turns out that from all the possible scenarios, the 
ones chosen by the best estimated values are precisely the two scenarios (Co > 0, Ci < 0) 
with Co + Ci < 3, and (0 < Co < 3, Ci = 0). 

The predicted values of the age of the Universe by the model when both (Co, Ci) are 
the free parameters, are a little smaller (11.72 and 10.03 Gyr) than the mean value of the 
observational constraint coming from the ages of the oldest galactic globular clusters but 
still inside of the confidence interval of this constraint (12.9 ± 2.9 Gyr). However, when 
we set Ci = 0, the age of the Universe evaluate at the best estimate of Co is 14.82 Gyr. 

We find that the computed minimum values of the x 2 -function by degrees of freedom 
(i.e., Xd.o.f.) are ver y dose to the value one (or even smaller to one), indicating that the 
"goodness-of-fit" of the model to the SNe la data is excellent. 

One of the disadvantages of the model parameterized by C = Co + (iH, is that 
the best estimated total bulk viscosity function ((z) is positive for redshifts z < 1 but 
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negative for z > 1. According to the local entropy production, negative values of £ 
imply a violation to the local second law of the thermodynamics. So, this is one of 
the drawbacks of the model. We find that the cause of this violation is due the fact of 
having considered the existence of the term (iH in the total bulk viscosity, this term 
characterizes a bulk viscosity proportional to the expansion ratio of the Universe. 

For the case where we set Ci = 0, we find that the total bulk viscosity is now 
positive and constant along the whole history of the Universe. This very simple model, 
in addition to the fact of that it does not violate the local entropy law, it has a good fit 
to the SNe la observations (measured through the Xd. .f.)> the age of the Universe is in 
perfect agreement with the observational constraints and it predicts a Big-Bang with a 
decelerated expansion followed with a transition to an accelerated epoch. This simple 
model of one parameter ( = ( is a better candidate to explain the present accelerated 
expansion than the model with two parameters C = Co + Ci-^- 

Finally, it is important to mention that to have viable bulk viscous models, a 
physical mechanism to explain the origin of this bulk viscosity must be proposed. There 
are already some proposals in that sense, for instance, a decay of dark matter particles 
into relativistic product is claimed as explanation of its origin. A preliminary study of 
this propose, using the SNe la test, have been done in [MJ [35] showing that it could be 
a promising mechanism to explain this origin. 
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